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Abstract 

We calculate the massive flavor non-singlet Wilson coefficient for the heavy flavor contri¬ 
butions to the polarized structure function gi{x,Q‘^) in the asymptotic region » m? 
to 3-loop order in Quantum Chromodynamics at general values of the Mellin variable N 
and the momentum fraction x, and derive heavy flavor corrections to the Bjorken sum-rule. 
Numerical results are presented for the charm quark contribution. Results on the structure 
function g 2 {x,Q‘^) in the twist-2 approximation are also given. 



1 Introduction 


Massless and massive contributions to the unpolarized and polarized structure functions in deep- 
inelastic scattering exhibit different scaling violations. For a precise determination of the QCD 
scale Aqcd or the strong coupling constant as(M|) their precise knowledge is therefore of im¬ 
portance [^. In the case of the polarized structure function gi^x^Q"^) the complete heavy flavor 


corrections are only available at 1-loop order [2,^ At higher orders in the coupling constant, 


the heavy flavor contributions were calculated in the asymptotic region based on the 


factorization derived in Ref. [^. Here denotes the virtuality of the exchanged gauge boson 
and m the heavy quark mass. The 0{a‘^) corrections in the polarized case were calculated in 
Refs. [^[^. In the case of the structure function the 1-Ioop heavy flavor corrections 

have been accounted for at next-to-leading order (NLO) QCD analysis [^. The corresponding 
flavor non-singlet corrections in the unpolarized case were calculated for pure photon exchange 
to 0{a1) in and in Ref. to 0(0^). 

In the present paper we calculate the O(a^) massive flavor non-singlet Wilson coefficient for 
the inclusive structure function gi{x,Q^) in the asymptotic region and also present 

the corresponding O(a^) result, extending Refs. (^[^. 

is given by 


The differential cross section for polarized deep-inelastic scattering 11-13 


(PaB 27ia^ 
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[Sl{x,y)gi{x,Q‘^) + Sl{x,y)g2{x,Q‘^)\ , 
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( 1 . 2 ) 


Here a = e2/(47r) denotes the fine structure constant, M is the nucleon mass, S = {p + /)^ is 
the center of mass energy of the lepton-nucleon system, with p and I the nucleon and lepton 
4-momenta, respectively, q = I — I' is the 4-momentum transfer and = —q^. x = Q‘^/{2p.q) 
and y = p.q/p.l are the Bjorken variables. denotes the degree of the nucleon polarization. 
The spin 4-vectors in the longitudinal and transverse cases are given by 


Sl = M(0,0,0;1) (1.3) 

St = M(0, cos(/3), sin(/3); 0) , (1.4) 


and ip denotes the angle between the vectors of the spin and the outgoing lepton. It contributes 
in a non-trivial way in the case of transverse polarization. 

The polarized structure functions are denoted by gpx^Q"^) and g 2 {x,Q^). In the leading 
twist approximation, the heavy flavor contributions to the structure function gpXjQ"^) is given 


by, cf. 14 


9i{x,Q^) 



Nf + 1 , 





^For an implementation in Mellin space, see 




[Afk{x, iP, Nf) + AfPx, iP, A^f)] 
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m^ jjL^ 
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(1.5) 


with the Np light flavor polarized (anti)quark densities, AG and AS = '^^J'^[Afk + Af^] 

the polarized glnon and singlet distribntions, and e* and eg the electric charges of the light 
qnarks and the heavy qnark Q, respectively, p denotes the factorization scale and ® the Mellin 
convolntion 


A{x) B{x) = / dxi / dx 2 S{x — XiX 2 )A{xi)B{x 2 ) ■ 

Jo Jo 

The actnal flavor non-singlet distribntion is defined by 

Nf 




k=l 


A/fc(x, pS Np) + A/fc(a:, pS Np) - — AS(x, pS Np) 

I\p 


( 1 . 6 ) 


(1.7) 


However, according to the representation (1.5), we will consider its whole first term, depending 


on as the non-singlet contribntion in what follows. The strnctnre fnnction g 2 {x,Q‘^) can be 


obtained from gi{x,Q‘^) nsing the Wandznra-Wilson relation 

The paper is organized as follows. In Section]^ we calcnlate the heavy flavor contribntions to 
the non-singlet Wilson coefficient in the asymptotic region to the strnctnre fnnction 

gi{x, Q"^) to 3-loop order in the strong conpling constant. We present the resnlts both in Mellin 
N and x-space. Nnmerical resnlts are given in Sectionj^ Conseqnences for the polarized Bjorken 
snm rule are discussed in Section and Section contains the conclusions. 


2 The Wilson Coefficient 


m 


The heavy flavor non-singlet Wilson coefficient contributing to the structure function gi{x, 
the asymptotic region 3> mJ receives its first contributions at 0{al). In previous analyses 
the tagged flavor case at 0{al) has been considered. In what follows we will refer to the inclusive 
case, i.e. the complete contribution to the structure function gi{x,Q‘^), and consider the terms 
due a single heavy quark. 

The non-singlet heavy flavor Wilson coefficient contributing to the structure function 
gi{x, in the asymptotic region S> mJ is given by 16 


+1) 




+ Os 


4 (2),NS 
^qq,Q 


(Nx + 1) + 


( 2 . 1 ) 


+1)++i)d,‘r'’(ivF+1)+ci?i'’‘"(ivF) 
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Here is the massive non-singlet operator matrix element (OME) and the label ‘Np + 1’ 

symbolically denotes that the OME is calculated at Np massless and one massive flavor, = 
as/( 47 r) = 5 f^/( 47 r)^ parameterizes the strong coupling constant, and we use the convention 

f{Np) = f{Np + 1) - f{Np) . (2.2) 


The calculation of the different contributions to the Wilson coefficient is performed m D = A + e 
dimensions to regulate the Feynman integrals. In the present polarized case the treatment of 
75 has to be considered. In the flavor non-singlet case both for the massive OMEs and the 
massless Wilson coefficients 75 always appears in traces along one massless line and there is a 
Ward-Takahashi identity which implies the use of anti-commuting 75 . 

The inclusive massive OME to 3-loop order for even and odd moments N has been 

calculated in Ref. 10 . The corresponding diagrams have been reduced using integration-by- 

H 


parts relations 17 applying an extension of the package Reduze 2 


The master integrals 

have been calculated using hypergeometric, Mellin-Barnes and differential equation techniques, 
mapping them to recurrences, which have been solved by modern summation technologies using 


extensively the packages Sigma 21,22 
and HarmonicSums (2^. 


EvaluateMultiSums, SumProduction 


psum 


The massless Wilson coefficients Og,g^(a:, from 1- to 3-loop order were calculated in 


setting the date terms in Cq^gf^^{Np) to zero, cf. also 


2- and 3-loop order were calculated in 5l 9 and 10 , respectively. 




Refs. 26|-|29 . At 3-loop order those of the structure f unction gi are obtained by that of F 3 29 

The non-singlet OMEs 


at 


For comparison, the massless flavor non-singlet Wilson coefficient in Mellin space is given 


by 28,29 


'1,NS 


(]Vf) = l + y]aJC<‘>7(iv^) . 


( 2 . 3 ) 


k=l 


In Mellin N space the Wilson coefficient can be expressed by nested harmonic sums Ss{N) 
which are defined by 

SbAN) = ^9 = 1, b, a,eZ,b,a,^ 0 ,N> 0 ,Ne N. (2.4) 

k=l 


In the following, we drop the argument N of the harmonic sums and use the short-hand notation 
Ss{N) = Ss- The Wilson coefficients depend on the logarithms 






and Lm = In ( ^ 


( 2 . 6 ) 


where the renormalization scale has been set equal to the factorization scale fi = = fip. 

As a short-hand notation we define the leading order splitting function Aygg^ up to its color 
factor 

3N^ + 3N + 2^ 


A7<? = 4 


2^1 - 


( 2 . 6 ) 


2N{N + 1) _ ■ 

The massive Wilson coefficient for the structure function gi{x,Q‘^) in the asymptotic region in 
Mellin space in the on-shell scheme is given by 


'h,NS 


(N) = i[i-(-i)q 



+ 41^7*? + i 
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'^The package Reduze 2 uses the packages Fermat 19 and Ginac 20 . 
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Here the color factors are given by Ca = Nc,Cf = {N^ — 1)/{2Nc),Tf = 1/2 in SU{Nc), and 

NS f3) 

iVc = 3 in the case of Quantum Chromodynamics. Cq^g^ (Np) denotes the massless Wilson 


coefficient at 3-loop order, cf. (2.2), and the polynomials Pj are given by 
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561iV^ + 1122iV^ + 767N‘^ + 302iV + 48 (2.20) 

1131iV^ + 2118iV3 + 1307Ar2 + 32N - 276 (2.21) 

1139iV^ + 2710iV3 + 635A^2 + 216iV + 828 (2.22) 

1199iV^ + 2398iV^ + IISIN^ + 18iV + 90 (2.23) 

1220iV^ + 2359iV3 + 1934Ar2 ^ 357JY _ ^33 ('2.24) 

3A^® + lliV^ + 10A^^ + 3A^2 + 7iV + 8 (2.25) 

12iV® + 16iV^ + 18N^ - - 5N - 8 (2.26) 

27N^ + 863iV^ + 1573iV3 + 1151iV2 + 144iV - 36 (2.27) 

51iV® + 102Ar^ + 121iV^ + 118A^^ + 48iV + 48 (2.28) 

648iV^ - 2103iV^ - 4278Ar^ - 3505Ar^ - 682N - 432 (2.29) 

1407iV'^ + 2418iV^ + 1793Ar^ + 13AN^ - 384Ar + 144 (2.30) 

-11145A^® - 32355iV® - 37b23N^ - 14329iV^ + 2392N‘^ + 120A^ - 1512 (2.31) 

-151N^ - A69N^ - 181Ar^ + 305iV^ + 208N^ + 40iV + 8 (2.32) 

3N^ + 9N^ + 70N^ + 77N^ + 39N^ - lOiV - 12 (2.33) 

6A^® + 18iV® - _ 20iV3 + AQN^ + 29N - 6 (2.34) 
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P28 = 15N^ + 36N^ + 30N^ + 8N^ + 3N^ + 16N + 20 (2.35) 

P 29 = 155iV® + 465A^® + 465iV^ + 371iV^ + 108A^^ + lOSiV + 54 (2.36) 

P 30 = 216iV® + 567A^® + 687iV^ + 381iV3 + 37iV2 - 44iV + 12 (2.37) 

P 31 = 309iV® + 807Ar® + 693iV^ - 271iV^ - 638iV2 + 68 Ar + 216 (2.38) 

P 32 = 525iV® + 1575iV® + 1535A^^ + 973iV^ + 536iV2 + 48iV - 72 (2.39) 

P 33 = 609iV® + 1485iV® + 1393A^^ + 83iV3 - 422iV2 + 156A^ + 216 (2.40) 

P 34 = 795iV® + 2043iV® + 2075Ar^ + 517iV^ - 298iV2 + 156iV + 216 (2.41) 

P 35 = 868 iV® + 2469iV® + 2487A^^ + 940iV^ + 27A^2 + 63iV + 72 (2.42) 

P 36 = 1770iV® + 4671iV® + 4765A^^ + 1205A^3 - 227A^2 + 1044Ar + 756 (2.43) 

P 37 = 7531iV® + 23673iV'^ + 23055iV^ + 7375iV^ + 1614Ar2 + 936Ar - 324 (2.44) 

Pgg = -4785Ar^ - 14355iV® - 4399iV'^ + 10327iV^ + 3548iV3 + SOOON^ 

+ 1080iV - 1728 (2.45) 

P 39 = 25iV^ + 138Ar6 + 311iV^ + 464Ar^ + 672iV3 + 670iV^ + 264iV + 48 (2.46) 

P 40 = -45iV® - 162iV^ - 858iV® - 1960Ar^ - 1885Ar^ - 1094iV^ - 804iV2 

-40iV + 192 (2.47) 

P 41 = 39iV^ + 138Ar^ + 847iV® + 1371iV® + 1283Ar^ + 485iV^ + lOliV^ 

+ 132iV + 72 (2.48) 

P 42 = 3549iV® + 14196iV^ + 23870iV® + 25380iV'^ + 15165iV^ + 1712iV^ - 2016iV2 

+144iV + 432 (2.49) 

P 43 = 5487iV® + 21948iV^ + 36370iV® + 28836iV'^ + 11943iV^ + 4312iV3 + 2016Ar2 

-144iV - 432 (2.50) 

P 44 = 10807Ar® + 43228Ar^ + 62898Ar® + 39178Ar® + 7027iV^ + 702Ar^ + 3240Ar^ 

+3456iV + 1620 (2.51) 

P 45 = 42591A^^ + 166764Ar^ + 245088A^® + 128254iV® - 26735A^^ - 40762Ar^ 

-3928A^2 - 1272Ar - 2160 (2.52) 

P 46 = -18351iV^° - 89784iV^ - 208773iV® - 267222iV^ - 192265Ar® - 46700iV® 

+ 14565A^^ + 7730iV^ + 1240iV2 + 1464Ar + 144 (2.53) 

P 47 = 165iV^° + 825A^® + 106856A^® + 321746iV^ + 396657iV® + 247433iV'^ 

+ 126914Ar^ + 51804Ar3 + 6336Ar^ + 4752iV + 5184 (2.54) 

P 48 = 828iV^^ + 7632Ar^° + 29217A^® + 59592Ar® + 66844Ar^ + 35738A^® + 7405iV'^ 

+ 16688A^^ + 27880A^^ + 11552Ar2 _ 33i2iV - 2304 (2.55) 

P 49 = 8274iV^^ + 78519iV^° + 313841iV^ + 686295iV® + 881001iV^ + 638778Ar® 

+204948A^® + 7992Ar^ + 32296A^^ + 26544Ar2 - 10656A^ - 8640 . (2.56) 


but 


We would like to note that we disagree with the 0(a^ ln(Q^//i^)) terms given in p 
agree with the representation in [^[5l] . 

One obtains the analytic continuation of the harmonic sums to complex values of N by per¬ 
forming their asymptotic expansion analytically, cf. 33,34 ^ Furthermore, the nested harmonic 
sums obey the shift relations 


SbAN) 


s,4n - 1 ) + 


sign(6)^ 

ivbi 


Sa{N) , 


(2.57) 


^These expansions can now be obtained automatically using the package HarmonicSums 
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through which any regular point in the complex plane can be reached using the analytic asymp¬ 
totic representation as input. The poles of the nested harmonic sums S's(iV) are located at the 
non-positive integers. In data analyses, one may thus encode the QCD evolution 35 together 
with the Wilson coefficient for complex values of N analytically and hnally perform one numerical 
contour integral around the singularities of the problem]^ 

over 


In x-space the Wilson coefficient is represented in terms of harmonic polylogarithms 37 


the alphabet {/o, /i, /_i}, which were again reduced applying the shuffle relations 38 . They are 
defined by 


Hb,ai^) = j dyfb{y)Hs{y), Ho^{x) =770 = 1, 

k 

fo{x) = -, fl{x) = f-l{x) = 

X 1 — X 


1 -I- X 


(2.58) 

(2.59) 


The Wilson coefficient is represented by three contributions, the (...)+-fnnction term, the 5(1—x)- 
term, and the regular term. Here the -I— distribution is defined by 


dy [F{y)]+ g{y) = / dyF{y) [g{y) - ^(1)] . 


(2.60) 


One obtains 

^NS 


< 1,91 


X = a. 


1 — X 


-CpTp 


2 + d^M 


'80 16 ■ 


-1- Ho 

[9 3 i 

+ Lq 


116 32 ^^ 

9 3 “ 


3 

3 


718 268 


27 9 


+ 5(1 - x) CfTf 


116 16 


— + —Ho + 8i7o' + — + -Ho 77i + -H^ + -Ho, 


8 . . 16 


9 


„rr 2 .21 .2 ^ 38 265 

2 [Lm + ~ ^ 


-\-C fTf 
+Lq 


--(x l) [Lq L 


M 


^(llx-l)-|(x + l)77o 


9 


^(I4x 5) y (x -h l)77o + ^{x + l)Hi 


__(218x + 47) 


■^(28x -7 13)i7o - 4(x + 1)772 + (^-^(I4x + 5) - ^(x + 


-^(x + 1)772 _ ^ + 1)^2 


+ai 


(1 -x)' 


;CaCfTf 


4 32 

-y (800x - 773)77 o 2 + -(94x - 12 l)C 2 


-hy(x-h2)i7o,i 

+^ 772 -^c; 

W 3 0 3 Si 


+ 


1 — X 


CaCfTf 


. 176 . 352 _ 

-Li—-Lj—+ L2 


27 


27 


M 


184 


+ L 


Q 


3104 704 16 , . 352 32, 

27 9 3 ° 9 3^ 


^For precise numerical implementations of the analytic continuation of harmonic sums see 
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+Lm 


1240 1792 248 ,,0 32 . 

81 27 9 ° 9 ° 


64 

IQHlH, + ?, 2 HoHo,i - -Ho,0,1 


)C2 + 96C3 


+ Lq 


80 . 30124 

- —Hn — 


9 


81 


14144 ^ 1216 

27 ° 9 ° 


27 9 3 °' 


4^ + yi7o)i^i -64i7oi^o,-i 


+ y 77 o - Vo,i - ^^^ 0 , 0,1 + 128 i 7 o,o,-i + 64 i 7 o,i,i 


192 + —i7o + 64i7ijC2-^C3 


43228 3256 496 , . 16 , . 


32 32 160 2 112 A 8 2 64 3 /368 

Y~Y “ 9 “ ° “ ^ ° j ' + 9 ° ^ ° ' ~ ° 


16 ^.2 , f o 128 ^ \ ^ , 128 A ^ 32 . 


+ -S-—Ho H, + —Hi )Ho,, - -Hi, + 


9 


1072 32 

-^- Ho 

27 9 ° 


320 


Hi \ Ho, 0,1 H—^ [T/q, 1 , 1,1 — -f 7 o, 0 , 1 , 1 ] + [ ——Ho — 32 Hi + 24 'j i 7 o,i,i 


160 


9 




496 


27 


-Ho 


112 


9 


160 


—Hl+ 8-—Ho]H: 


9 


9 1^9 0,1 IS2 ^ 1 27 


9 9/3 3 




CM 

Ot 

+ CITf 


16 - Ho - Hi 

[ 3 ° 3 J 


+LmLq 
16 


1 ,-^^Ho-^-^Hi 

3 ° 3 


64 


— Hi + 1677 / - —C 2 
3 ° ^ 3 ^ 


+ L\ 


+ 

334 


Q 


( 128 
8 H ——Ho 1 H 


9 


9 


+ LmLq 


88 176 

- Ho + 

3 9 ° 


9 3 ° ' 


3 0 3 Si 


64 


^-^Hl + ( ^ + 


152 


64 

- —77i )/2 - 64/3 


+ Lm 
1424 


206 112 


-^Ho + -Hl+\— + -Ho]Hi 


9 

+ Lq 


3 3 

160 

Lo 
3 ° 

2360 4508 


160 32 


1 64 

77o + ^77nM i/i - + 


3 
784 


9 3 ° 


9 


^ 27 3 9 9 


4556 3680 


27 


9 


Ho-l28Hl ]Hi + 


512 


— 12877o 1 H^ + 128HoHo,—i 


64 


64 


48 — 3 “-^o + “^-^1 )-^o,i — 25677o,o,-i — 6477o,i,i + 


/2608 832 


448 


3 
302 

~V 

128 


Hi 1/2 + 320/3 
13624 


14197 

54 




27 


27 


81 


81 


-Ho 


1628 ,,2 304 

-Hi + —1 


27 


27 


HoHi+[—- 


112 1304 


27 

32 


^^0 Vi + f ^ - ^^^0 V ? 


V 9 


160 


9 

128 


—Ho - -jHl + -^HoHi - —Hi ]Ho,i 
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16 2 128 
9 9 


64 


0,0,0,1 + -^-f^0,0,l,l + 


2488 1192 


3 J 9 ^ 9 ’/ 


664 


27 

\ , /3088 

. 1 K 2 + 


27 


TT Ij 2 I 

^0 “ + 


160 

~9~ 


V 27 9 9 3 


464 



o 128 o 64 , 

'320 64 

- 

+ CfT^ 

Ll^ + Ll- + Ll, 

W + 

+ Lq 


9^9 


1984 

+ 


3760 2144 

+ 


81 


27 


27 


928 128 


^ 0 + — + —Ho\Hi 


3 ° 9 ^ 9 9 ^ 


27 

160 


9 


12064 64 160 . 32 ^ „ 

729 81 81 ° 81 ° 


896 

+Lm 

256 


Cs 


+ CpNFn 




27 


2176 


320 32 

■-On- —■tin 


27 


9 


Q 27 ^ 

■7520 


+ Lq 


81 


H \H +^28„2 , 256 512 


+ 


+ 


928 

4288 


9 ° 9 ^ 


128 2 , A856 
—Ho + —ifo + 

24064 128 320 ,,2 

i7n - —i7n 


_64 3 _^ 
81 27 



-fC3' 


+ 


938 16 


+ 5(1 - x) I CaCfTf 
1595 


729 ^ 81 
44 


81 


~^M~F, - ^Q~Fr + ^'m 


9 


32 1024 196/ 


+ Lm 
55 


272 68 

27 +^‘'=’ + y‘'^ 


9 

+ Lq 


34 

y 

11032 

27 


+ 


243 


10045 16 2624 176 

81 ^3 “ ■/ 82 C 3 + - —Cs 


-8 Ba 


+ CITf 


-\-Lm Lq2 + L 


M 


-5 


-10 + —Cs 
3 ^ 


+ 


-48 + —Cs 
3 ^ 


112 136 

-^Cs - ^C4 


+ Lq 


368 64 

y + y^^ 


1616 

~9~ 


Cs 


392 


C 4 



2039 

yy 

32 


13682^ 32^^ 3304^ 352^ 

^^Cs + yC 2 Cs - ^C 4 + 1654 + —Cs 


16 


152 


496 


y 77 + -^QTT + - Lq— + Lai^ + Lq 


1624 3658 


g . -Q g Mg Q" g ■ 27 ' ^ 27 


243 


16 


32 


304 


700 


y 7” + ~ ^QHT + 


27 


3248 4732 

+ 


27 


243 


176 


+ l) + + l) + L 


--(83X- 37) 


+f (x +1]H„+ [26 - Hi] 

3^ 3x + l'- 3x + l'- 
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+-^Q 

176 

iT 


_1 (865a: + 109) - (x + l)ff, + - H_,H„] 


X 


l)H,+ 


16 X 
3 a; + 1 


[ 2(2 - HS] 


+ Lm 


■ — (4577a;-4267) 
81 ^ 


16 2 X 419a;2 + 4a; + 25 , . 32 a; 

—-(29a; - 109) i7o +- Hi - 

27^ ^^9 a; + l ° 9a; + l 


+8(a; + l)i7g^ 


Hi + 


128 4a;^ + 3a; + 4 


9 


+ 


128 4a;^ + 3a; + 4 


9 a; + 1 
32 a;2 + 1 


Hn 


a; + 1 
16a;2 + l 


+ ( -y 


- 16//o) {x + l)Hf 


3 a; + 1 
64 a;2 + 1 
3 a; + 1 
" 4 


3 a; + 1 

[Hofl^-i — 2i7o,_i,i — 2i7o,i,_i] + 

16 3a;^ + 4a; + 3 


[4i7o,i-^o] )i7-i + -—77. 


77-1 


a; + 1 


Hq C 2 — 32 


3 a; + 1 

/16 3a;2 + 14a;-9 
\ 9 a; + 1 
a;^ + 3a; + 1 


0 , 0,1 


a; + 1 


Cs 


+7q 


(12329a; - 577) + f 


f 64 181a;2 + 239a; + 49 32 (a; - l) ^ 


27 


a; + 1 


32 6a;^ + 25a;3 + 18a;2 + 25a; + 6 
+ —--77-1 77 o + 


9 


X 


+ 1 ) 


X 


H\ 

3 a; + l 

8 12a;3 - 21a;2 - 77x - 24 

9 a; -f 1 


16 5a;2 - 2a; + 5 


+ 


3 

352 


a; + 1 


77-1 ]Hl + 


80 X 
9 a; + 1’ 


773 + (-(703a;+ 253)--(a;-3)77, 


9 


(a; + 1)77o^ 77i + (^^ (a; + 1) - y (a; + 1)77o^ 772 


+ 


/64 3a; + 1 
It a; + 1 


77n 


32 6a;^ + 25a;3 + 18a;2 + 25a; + 6 64 (a; - l) ^ 


(a; + 1 ) 


X 


3 a; + 1 


77_i 77o,_i+ — (a; + l) 


208 


16 


64 



32 


+ —(a; — 3)77o + —{x + l)77i |77oy + —{x + 3)77o,o,i ~ 32(a; + l)77oyy 


64 (a; - l)' 


77o,-i,-i 


32 5a;2 + 10a; + 9 


3 a; + l ’ ’ 3 a; +1 

16 12a;3 - 23a;2 - 72x - 17 32 (a; - l) 


/ 64 

77o,o,-i + ( ~“^(^ + 2)77 i 


9 a; + 1 

64 3a;^ + 3a; + 2 
'Y a; +1 
32 199a;2 + 174a; + 199 


3 a; + 1 


-77_i- 


16 3a;^ + 8a; + 3 
3 a; + 1 


77o C 2 


2 / 4 

— (108295a; - 8668 I) + - — (995a; - 2807) 
729 ^ y 81^ 


R I I 64 a;^ + l 3 . 


81 a; + 1 

4 253a;2 + 391a; + 586 16 19a;2 + 18a; + 19 


81 a; + l 27 

8 22a;2 + 7a; + 25 32 a;^ + 1 


81 


a; + 1 


27 a; + 1 


a; + 1 

H-iVo- 


56 


77-1 + 


16 a ;2 + 1 

9 a; + 1 


77!i 77, 


16 X 
Ya; + 1 


H« + (^{9x + 4)H„ 


^^(65x-29) + l(Ux + 3)H-^+-^^{x + 


27 


27 



■-(2a; + 5)77o--(a; + 



77i + (^--(43a; - 46) 
77i + ^(a; + l)77o77f + + l)77_ 


14 



32 199x2 + 174a: + 199 , 64x2 + 1 3 , /^2564x2 + 3x + 4 


81 


X + 1 


9 X + 1 


Hi, ]Ho,-,+ 


{ 


27 


X + 1 


H. 


-1 


— — (l43x + 2 ) —^ (l3x + 6) Ho — -[x + 1 ) Hq + (llx + 20) 

+ “^(^ + ^)Ho^Hi —^(x + l)i72^i7o,i + 


128x2 + 1 


77_i i7o,-i.-i- 


256 4x2 + 3x + 4 


27 


9 X + 1 

32 19x2 + 18x + 19\ 

- j i7o,o,-i + 


X + 1 




64 x2 + 1 
9 X + 1 


H-i 


27 

16 


X + 1 


8 9x2 + 101x + 12 64x2 + 1 ^^ 

+ —— —H_i 


-^{7x + l)Ho-^{x + l)H,jHo,o, - 27 


+ 


(-y(^ + 4 


160 

IT 

128x2 + 1 
9 X + 1 


27 x + 1 '9x + l 

256 4x2 + 3x + 4 


x + 1 


-Ho, 1,-1 


16 


H-i —— (llx + 5)i7o + 16(x + -^ 0 , 1,1 



64 x2 + 1 
9 x + 1 


+ 2-f7o,- 1 , 1,1 + Hofl ,- 1,-1 ~ -f^o,o,-i,i + -f^o,o,o,-i 


~-f^o,o,i,-i + 2i7o,i,-i,i + 2-f7o,i,i,-i] + 


64 5x2 ^Qx-i 


9 


H 


0 , 0 , 0,1 


16x2-2x-11 ^ 112, 

-ffoAl.l - -g-(a:+ + 

32 29x2 + 18x + 29 
x + 1 

8 63 x 2 _|_ 292; 5 


x + 1 

/16 174x2 + 209x- 189 


V81 


x + 1 


27 


x + 1 


H-i + -(3x + 14 ) + — (x + l)i7o^-f^i 

32x2 + 1 ^ \ ^ 32x2 + 1 „ 64, ^ 

-i7-i i7o- H\ H-(x + 1)77 

9x + l 7 9x + l"^ 9^ ^ 


8 3 x 2 _j_ 32 : + 9 


+ 9 


x + 1 
128x2 + 1 


772 


— (7x + 1 ) 77 o 4 ^ C 2 + ^ 


2 497x2 + 1102x + 1085 


27 


x + 1 


77_i 


32 6x2 + 42 ; _ 3 


9 X + 1 9 

8 36x2 ^ ^ 22 


x + 1 


C 4 


x + 1 


+ CITf 


16 \ 16 
77o + — (x + 1)77 iJC 3 + — (x + 1)774 


[^Q + ^m^q] 


32 


--(x + 5) + —(x + l)77i 


+8(x + l)77o 


+ L 


M 


28(2x- 1) + (^-^(llx + 5) + 


3^ '3 

64 x2 + 1 


77-1 77, 


4 9x2 + lOx + 9 ^^2 / 16 64, 

3 x + 1 ° V 3 ^ ^ 3 ^ 

i\rr 2 64x2 + 1 8 9x2 + lOx + 9 

+ HI - T7TT"»-‘ + 3 . + 1 



3 x + 1 


0 I -f^i ~ 2 (^ + l) Hq^i 


+ 


-(I61x + 130) 
9 


+ |^-y (l5x + 4) + y ^ )^^^ 77-i^77o - 24(x + 1)772 |^_^^ 5 g 2 . 17 ) 


160 


X 


3 

64 x2 + 1 
3 X + 1 



4 21x2 + 34x + 21 ^ „ 8 23x2 + 38x + 23 , 

0 1 77 i- 77 n H-C2 

3 x+1 °3 x+1 ^ 


77o,-i — 8(x + l)77o,i 


+ LmLq 


- (I9x - 85) + - (l3x + 1 )Ho 
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4 195x2 + 238x + 123 ^ 32 3^^ + 4x + 3 ^3 


X + 1 


+8{x + l)H^+ |^^(4x + l) + y(x +l)i/o^^?'i - y(a^ + l)C 2 

, r r 4 4195x2 + 238x + 123 , 32 3x2 + 4x + 3 . 

+Lm^ g(337x + 235)ffo g H, ^ x + 1 ^ 

+ (^-^ (287x - 113) - ^(5x + 2)i7o - y (x + 1) H, + ^{l78x- 125) 

+ (^-y(7x + 3) - y(x + l)i7o^^i + (^^(x + l) + y(x + l)i7o^77o,i 


256 4x2 + 3x + 4 
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16 3x2 _ 2a; + 3 


X + 1 


256 4x2 + 3x + 4 
9 X + 1 


-fixT + y [2"o.-u - -ffoA-i + 2i/o,.-i] 

/8 117x2 + 118x +81 128x2 + 1 16 (x + 3) (3x + l) 

\9 x + 1 3 x + 1 3 x + 1 ^ 
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y x+ 1 
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-- (12332X - 4905) + - (a; + l)i7!, + - 


+ —(- 10999X- 8399) + — 


128 x2 + 1 


27 x + 1 


32 21x2 + 30x + 17 
~ y x + 1 
64 199x2 + 174x + 199 ^ 


x + 1 

32 19x2 + 18x + 19 


^-1 ^0+ Uz 


32 x2 + 1 2 2 4179x2 + 5255x + 2868 


- —H\ - 

9 x + 1 81 

_ 10 177x2 + 218x + 105 
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,,, 1 51x2 + 70x + 51 
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° 27 x + 1 
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+ + (-^(+ + 1 )+ 
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(2.61) 


Again, we used the short hand notation Hs{x) = 77^ also here. The transformation of the Wilson 
coefficient to the MS scheme for the heavy quark mass affects the massive OME at 3-loops and 
was given in Ref. 10 ; the terms are the same in the unpolarized and polarized case. 

The non-singlet contributions to the structure function g 2 {x,Q'^) can be obtained via the 


Wandzura-Wilczek relation 15 


92 {x,Q^) = -gi{x,Q^) + [ —gi{y,Q^), 

Jx y 


(2.62) 


where both structure functions refer to the twist-2 contributions. This relation is implied by a 
relation of the OMEs in the light-cone expansion, cf. 
in the covariant parton model in Refs. 


39 . The relation has also been proven 


40 42 . For gluonic initial states, it was derived in 43 


Eq. (2.62) also holds including target mass corrections 44,45 and hnite light quark contributions 


corrections 47,48 


45 . Furthermore, it holds in non-forward 46 and diffractive scattering, including target mass 


3 Numerical Results 


In what follows, we will choose the factorization and renormalization scale = Q^. We hrst 
study the behaviour of the massive and massless Wilson coefficients in the small and large x 
region and then give numerical illustrations in the whole x-region. 

At small X, the pure massive Wilson coefficient behaves like 


+ oc al4CpTp\n\x) + a 
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ln^(l — x) 
1 — X 


(3.2) 

There is a term oc ln^(l — x)/{l — x) at 0(ln(Q^//i^)), being of relevance for different choices of 
the factorization scale. 

The above results can be compared with the case of the massless Wilson coefficient 


10 


K al-CFTpln\x) 


(]Vf) oc a: 




ln^(a;) 
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oc al-CpTp 
CTifHNx) oc 


ln^(l — x) 

1 — X 
ln^(l — x) 
1 — X 


(3.3) 

(3.4) 

(3.5) 

(3.6) 


The small x behaviour can be compared with leading order predictions for the non-singlet evo¬ 
lution kernel in Refs. 49, ^ . Indeed both the massive and massless contributions follow the 

ln^^(a;). However, as is well known less singular terms widely 


CkOj 


principle pattern 

cancel the numerical effect of these leading terms. For the large x terms the massless terms 
exhibit a stronger soft singularity than the massive ones. 

In the following numerical illustrations we use the polarized parton distributions of Ref. [^, 
which are of next-to-leading order (NLO), since no next-to-next-to-leading order (NNLO) data 
analysis based on the anomalous dimensions calculated in Ref. has been performed yet. The 
values of correspond to those of the unpolarized NNLO analysis 52 . The heavy and light 


flavor Wilson coefficients being discussed in the following are given in Eqs. (2.2) and (2.3) 


In Figure 1, the 2- and 3-loop heavy flavor corrections to the non-singlet term of the structure 
function xgi^x^Q"^) are calculated in the case of charm, assuming rric = 1.59 GeV 53 , using 
the formula for the Wilson coefficient Eq. (2.61), and setting /i^ = Q^. With growing 


the distribution diminishes at larger values of x and grows towards medium values. The O(a^) 
corrections lead to stronger effects if compared to those at O(a^). We have applied the asymptotic 
Wilson coefficients for all the values given here, which only holds for values jm} > 10. For 
the heavy quark distributions we formally show also the result at = 4 GeV^, outside this 
region, indicated by dotted ( 0 ( 0 ^)) and dash-dotted lines ( 0 ( 0 ^)). 

Figure 2 shows the effect of the Wilson coefficients comparing the contributions from 0(a°) 
to 0{a^) at = 4 GeV^ as an example, where a depletion is obtained with growing order. The 
3-loop light flavor contributions to xgi{x,Q^) {Np = 3) are illustrated in Figure 3. Here the 
evolution is strengthened by growing in the large x region and depleted for lower values of 
Q^, considering only the effects due to the Wilson coefficient. 

In Figures 4 and 5 we illustrate the ratio of the flavor non-singlet charm corrections to those 
by the light quarks given in Eq. (2.3) up to 0{a^) and O^a^), respectively. At 0{a‘l) the effect is 


of 0(1%) and below, for the lower scales Q^, but higher values are obtained for very large scales 
as ~ 1000 GeV^ in the region x ~ 0.003. A qualitatively similar picture is obtained including 
the O(a^) corrections. The effect on the ratio is about doubled. To resolve relative 

effects of 0(2%) requires higher luminosities than available in present day experiments. They 


may become available in the planned experiments at a future EIG 54 


Figure 6 shows the 2- and 3-loop charm flavor non-singlet contributions to the structure 
function xg 2 {x,Q‘^) according to the Wandzura-Wilczek relation (2.62) implying the oscillatory 
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behaviour. In size these effects are comparable to those of the structure function xgi{x,Q‘^) 
shown in Figure 1. With growing the effects become somewhat smaller. In Figure 7 we show 
the corresponding massless contributions to the structure function g 2 {x, Q"^) at = 4 GeV^ for 
the different orders in a^, which slightly diminish adding higher order contributions. Taking into 
account the O(a^) corrections, the light flavor corrections to g 2 {x, Q'^) ( 1.5[2.62 ) grow somewhat 
in size with larger values of Q^, see Figure 8. Similar to the case of the structure function xgi 
the 0 {al) charm flavor non-singlet corrections to the structure function xg 2 {x,Q‘^) amount to 
0 ( 1 %). 


4 The Bjorken Sum Rule 


The polarized Bjorken sum rule 55 refers to the first moment of the flavor non-singlet combi¬ 
nation 


dx [c/f (x, Q^) - gl^{x, Q^)] = - 


9a 


9v 


OBj(as), 


(4.1) 


with gAy fhe neutron decay constants, gA/gv ~ —1.2767 ± 0.0016 56 and ds = otsj'x. The 


1- | 57 |, 2- 158 ], 3- 30 and 4-loop QCD corrections |31] in the massless case are given by 


C'Bj(a,), = 1 - a, + a2(-4.58333 + 0.33333Nf) + a3(-41.4399 + 7.60729A^j. - 0.17747iV|) 


+at(-479.448 + 123.472iVi. - 7.69747iV| + 0.10374Ar3' 


(4.2) 


choosing the renormalization scale fy = cf. 28 for SU{3)c- Here Np denotes the number of 


active light flavors. The expression for general color factors was given in Ref. 31 


For the asymptotic massive corrections (2.2) only the first moments of the massless Wilson 


coefficients Cq^gf^’^' ^{Np) contribute, since the first moments of the massive non-singlet OMEs 
vanish due to fermion number conservation, a property holding even at higher order. Therefore, 


any new heavy quark changes Eq. ( |4.2[ ) by a shift in Np -A- Np + 1 only, for the asymptotic 

at O(a^), where no 


corrections. Different results are obtained in the tagged flavor case 
inclusive structure functions are considered. Corresponding power corrections were derived in 


59,60 


5 Conclusions 


We calculated the heavy flavor non-singlet Wilson coefficients of the polarized inclusive structure 
function gi{x,Q‘^) to 0{al) in the asymptotic region ^ m^. The first contributions of this 
kind are of 0 (q!^). In the case of twist-2 operators the corresponding contributions to the 


15 


structure function g 2 {x,Q'^) can be obtained using the Wandznra-Wilczek relation (2.62) 
cf. 39 42,45 . The asymptotic Wilson coefficient is obtained by using the factorization formula 
[^, Eq. (2.2), based on the massive OME [To] and the massless Wilson coefficient 29 to 3-loop 


order. The heavy flavor Wilson coefficient can be thoroughly represented by nested harmonic 
sums in Mellin-77 space and by harmonic polylogarithms in x-space. We presented numerical 
results correspo ndin g to the charge weighted polarized parton contribntions oc A/(x, Q^) -|- 
A/(x, Q^), cf. (1.5), referring to the polarized parton distribution fnnctions at NLO for an 
illustration. Comparing with the corresponding massless cases the heavy flavor corrections in 
case of charm are of 0(1 — 2%), requiring high luminosity experiments to be resolved, which 
are planned for the future electron-ion collider EIC 54 . We also considered the contribution 
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of the asymptotic Wilson coefficient to the polarized Bjorken sum-rule. Due to fermion number 
conservation for the massive flavor non-singlet OME in all orders in ag, only the hrst moment of 
the massless Wilson coefficient contributes and the effect of each heavy flavor results in a shift 
of Np by one unit in the expression for the massless polarized Bjorken sum-rule. The results of 
the present calculation could be easily applied to derive the asymptotic heavy flavor corrections 
to the neutral current structure function xG^,-, |^. However, the corresponding massless Wilson 
coefficient to 3-loop order has not been calculated yet. 
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Figure 1 : The 2- and 3 -loop non-singlet charm contributions to the structure function xgi{x,Q‘^) by the 
asymptotic heavy flavor Wilson coefficients in the on-shell scheme for nic = 1-59 GeV. Here we used 
the value of q;s(M|) = 0.1132 and the NLO parton distribution as reference. Figures 2-8 below are 
calculated using the same setting. 



Figure 2: The light flavor contributions {Np = 3) to the non-singlet charm contributions to the structure 
function xgi{x,Q'^) at = 4 GeV^ illustrating the contributions for the different orders in Og. 
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Figure 3: The light flavor contributions {Np = 3) to the non-singlet charm contributions to the structure 
function xgi{x,Q'^) at 0{ag) for different values of Q^. 



Figure 4: The ratio in the non-singlet case at 0{al) for different values of Q^. 
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Figure 5: The ratio 


^charmy'^iightj^^ non-singlet case at 0{a^) for different values of Q^. 



Figure 6: The 2 - and 3-loop non-singlet charm contributions to the twist 2 contributions of the structure 
function xg 2 {x,Q‘^) by the asymptotic heavy flavor Wilson coefficients in the on-shell scheme. 
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Figure 7: The light flavor contributions {Np = 3) to the non-singlet charm contributions to the structure 
function xg 2 {x,Q'^) at = 4 GeV^ illustrating the contributions for the different orders in Og. 



Figure 8: The light flavor contributions {Np = 3) to the non-singlet charm contributions to the structure 
function xg 2 {x,Q^) at O(a^) for different values of Q^. 
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